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Abstract. We study the Ising model two-point diagonal correlation function 
C{N,N) by presenting an exponential and form factor expansion in an integral 
representation which differs from the known expansion of Wu, McCoy, Tracy 
and Barouch. We extend this expansion, weighting, by powers of a variable A, 
the j'-particle contributions, fj^ j^- The corresponding A extension of the two- 
point diagonal correlation function, C{N,N;X), is shown, for arbitrary A, to be 
a solution of the sigma form of the Painleve VI equation introduced by Jimbo 
and Miwa. Linear differential equations for the form factors ff^ j^f are obtained 
and shown to have both a "Russian doll" nesting, and a decomposition of the 
differential operators as a direct sum of operators equivalent to symmetric powers 

(7) 

of the differential operator of the elliptic integral E. Each fj^ pf is expressed 
polynomially in terms of the complete elliptic integrals E and K. The scaling 
limit of these differential operators breaks the direct sum structure but not the 
"Russian doll" structure. The previous A-extensions, C{N,N; A) are, for singled- 
out values A = cos{7rr7i/n) (m, n integers), also solutions of linear differential 
equations. These solutions of Painleve VI are actually algebraic functions, being 
associated with modular curves. 
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1. Introduction 

The two dimensional Ising model in zero magnetic field is, historically, the most 
important solvable model in all of theoretical physics. The free energy [1], the 
partition function on the finite lattice [2] and the spontaneous magnetization [31 H] 
were computed long ago by Onsager, Kaufman and Yang. These computations, and 
subsequent studies of the correlation functions [5]- [9], form the basis of scaling theory 
and of the renormalization group approach to critical phenomena. 

The next most important macroscopic property of the Ising model, which one 
would like to compute, is the magnetic susceptibility at zero magnetic field, which is 
expressed in terms of the two-point correlation functions C (Af, N) with a spin at the 
origin and the other spin in row M and column N , as 

oo 

kBT-x= E {C{M,N) ~M\Q)) (1) 

J\/,JV=-oo 

where A^(0) is the spontaneous magnetization (which is only non zero for T < Tc). 
Unlike the free energy, and spontaneous magnetization, this has no known closed form 
expression, and the study of the magnetic susceptibility has been the most challenging 
outstanding problem in the field for over 50 years. 

The first serious analytic study of the susceptibility was made in 1976 by Wu, 
McCoy, Tracy and Barouch [9j who used their expansions of the correlation functions 
to write the susceptibility as an infinite series in multiparticle contributions as 

oo 

X±{T) = E ^^'^ (2) 
j 

where the subscript ± refers to T above (resp. below) and the sum is over odd 
(resp. even) values oi j for T above (below) Tc. In ref. [S] the contributions x*"^"* a-nd 
-^^gj-g gxplicitly calculated. 

No further analysis of the susceptibility X; or of the x^"'^ j was attempted until 
1999 when Nickel, in two remarkable papers [10l[Tl], showed for j > 3 that the 
have singularities in the complex temperature plane whose number increases with j 
and become dense on a circle j oo. Unless a remarkable cancellation takes place 
this discovery implies that the magnetic susceptibility will have a natural boundary in 
the complex temperature plane which extends to Tc- This natural boundary is a new 
phenomenon which is not incorporated into scaling, or renormalization theory, and, 
thus, it raises significant questions about our understanding of critical phenomena. 
Consequently it is most important to deeply understand the properties of the Ising 
susceptibility, and this challenging question certainely requires some serious progress 
on the two-point correlation functions of the Ising model. Note that some re-summed 
high temperature series |121I13| in the anisotropic case has, already, enabled Guttmann 
and Enting |13j to conjecture, for the anistropic x, a natural boundary in one variable 
when the second variable is fixed. 

In 2001 the work of Orrick et al. [T3] provided a polynomial time algorithm 
for obtaining the coefficients of the susceptibility series of the two-dimensional Ising 
model: from a combinatorial enumerative viewpoint this can be viewed as a "solution" 
of the problem. The existence of such a polynomial time algorithm for a lattice 
problem, instead of the exponential growth of the calculations one expects at first 
sight, can be seen as some "combinatorial integrability" of the model [T3]. However a 
(very) efficient way of getting very large series expansions for a physical quantity of 
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a model of lattice statistical mechanics, is far from providing the closed formula and 
exact results one might desire : for instance, there is still a lot of work to be done in 
order to extract singular points, singular behaviours, from the knowledge of very large 
series. 

In 2004 several of the present authors [151 Ull IIZl HI] initiated the study of the 
Ising susceptibility, beyond the singularity analysis of Nickel [101 HI] i by determining 
the Fuchsian linear differential equations for x''^'' ^-^d x^^'' ^ ^ function of the 
temperature. These equations have many remarkable properties such as a "Russian 
doll" nesting structure : the function x^^^ satisfies the equation for x^"^^ and x^^^ 
satisfies the equation for x*"^"* • If this nesting can be proven to extend to all of the 
X*--'-' 's there must be remarkable structures in the Fuchsian equations and the hope is 
thus raised that it may be possible to characteriz^ the full susceptibility. 

In a more recent paper [19j several of the present authors provided new results 
on the exact expressions of the two-point correlation functions of the Ising model, 
especially the diagonal correlation C{N,N), underlining the key role played by the 
second order linear differential operator corresponding to the complete elliptic integral 
of the first or second kind K oi E. 

In this paper, we study the diagonal correlation functions C{N,N) as a form 
factors expansion. Our starting point will be the expansions of the diagonal 
correlations in an exponential form [9J, both for T < Tc 

oo 

C.{N,N) ^ (1 - . exp(^ 4^^'^) (3) 

n=l 

with 

t = (smh{2Ey keT) smh{2E'' /kBT)^ ^ (4) 

and for T > 

OO OO 

C+iN, N)^ (1 - t) V4 . J2 G^^:^'^ ■ exp [Y. ^i^+U+i) (5) 

n— n— 1 

with 

t^ ((sinh(2£;7fci3r)sinh(2£;'VfcBT))^ (6) 

where E^ and E"" are the horizontal and vertical interaction energies of the Ising 
model. When the exponentials in ^ and ([5]) are expanded, the correlations can also 
be written in what is called a "form factor" expansion : 

oo 

C.{N,N) = (1 - . (l + ^ /(.^S) (7) 

oo 

C^{N^N) = (l-O^/^-E/K"/'^ (8) 

71=0 

The form factor fl^\ is interpreted as the "j-particle" contribution to the two-point 
correlation function. It is natural to consider A-extensions [31 [3T] of the previous 

§ The full susceptibility could be the solution of a nonlinear equation, or the solution of a system of 
PDE's, or solution of a nonlinear functional equation, or ... 
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OO 



C_(iV,iV;A)= [l^tf/'-il + yX^-f^^'^il) (9) 



functions 

n (AT AT- \\ — l^ _ / 1 I , 

/ . ^ Jn,n 

n=l 

oc 

C+{N, N; A) = (1 - . ^ • (10) 

which weight each fj^^pj by some power of A, and to interpret A as being analogous to 
a coupHng constant in a quantum field theory expansion. Such A-extensions naturally 
emerge from the Fredholm determinant framework in Ql] . We will present new integral 

representations for F^^^, ^^n^n^'' ^-^^ /n^'n sec.©- We will see that they are 
much simpler, and more transparent, than the forms obtained from C{M,N) of ^ 
by specializing to M = iV. The proof of these results is obtained by extending the 
expansion solution for the leading term given in 1966 by Wu 7 , to all orders. It will 
be published elsewhere. 

The diagonal correlations C{N, N) have the property, discovered by Jimbo and 
Miwa [20] in 1980, that their log-derivatives are solutions of the "sigma" forni of a 
Painleve VI function 

,da f , .da 1\ f da \ , _ , 

-'Tt[^'-'^Tt-^-A)['Tt-') ^''^ 

where a is defined for T < Tc as 

, . dlnC-(N,N) t , , 

aN{t) - tit~l) -- (12) 

with the normalization condition 

C-{N,N) = 1 +0(t) for t^O (13) 
and, for T > Tc, as 

, , dliiC+(N,N) 1 
a^it) = t{t-l) ' -- (14) 

with the normalization condition 

C+{N,N) = 11/^ . . (1 + C'(i)) for t^O (15) 

where {a)N — r{a + N)/r{a) denotes the Pochhammcr symbol. 

One can easily verify that (fTT|) . the iV-dependent sigma form of Painleve VI, is 
actually covariant by the Kramers- Wannier duality : 

{t,a,a',a") (1, ^, a - 1 ■ a' , ■ a") (16) 

On another hand, Jimbo and Miwa introduced in 20j an isomonodromic X- 
extension of C{N, N) and showed that this more general function C{N, N; A) also 
satisfies (|lip. The motivation of introducing an isomonodromic parameter A, in 
the framework of isomonodromy deformations, is, at first sight, quite different from 
the "coupling constant" motivation at the origin of the form factor A-extensions © 
and (jlOp . In sec. ([3]) we show that these two A-extensions are actually the same by 

X We use a variable t which is the inverse of the one of Jimbo and Miwa | 20| . 
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demonstrating that the recursive solutions of ITTI) , analytic^ in i , agree with and 
(fTO]) where the /^'^^'s are obtained from C±{N, N; A), the A expansion of C±{N, N) 
of sec. (121). The normaHzation condition (fT3|) fixes one integration constant in the 
sohition to (fTTj) . We find that the second integration constant is a free parameter, 
and, denoting that parameter by A, that our one parameter family of solutions for 
(N, N) can be written in a form structurally similar to the right hand side of ([9]) . 
Furthermore, we have confirmed, by comparison with series expansions of the multiple 
integral formulas for /^\r derived in sec. ([2]), that this family of solutions is, in fact, 
identical to C-{N, N; X) as defined in ([9]). Similarly, the condition (fTS]) gives rise to 
a one parameter family of solutions for C+{N,N) that is identical to (fTO)) . After all, 
the fact that these two distinct A-extensions of C± {N, N) identify is not altogether 
surprising, since Jimbo and Miwa's derivation of (|lip also starts from a multiple- 
particle expansion of the correlation functions in terms of free fermion operators. It 
does not, however, appear to have been observed previously. 

In secdH) we use formal computer algebra to study the functions /^^^- We obtain 

the Fuchsian linear differential equations satisfied by the /^^^ for fixed j < 9 and 

arbitrary TV. We also find the truly remarkable result that the families /^'^f^^'' and 

fN'^N each annihilated by linear differential operators which have a nested "Russian 
doll" structure. Beyond this "Russian doll" structure, each linear differential operator 
is the direct sum of linear differential operators equivalenlQ to symmetric powers of 
the second order differential operator corresponding to /^^^, (or equivalently to the 
second order differential operator Le, corresponding to the complete elliptic integral 
E). A direct consequence is that the form factors f^^''^^\ and f^% are polynomials 
in the complete elliptic integrals of the first and second kinds, K and E: 

K = (1/2, 1/2; l-t), E ^ 2F1 (1/2, -1/2; 1; t) (17) 

A simple example is /q — K ■ {K — E)/2. 

The closed formula we obtain for the differential operators in these nested 
"Russian doll" structures, enable us to take the scaling limit of these operators. We 
study this scaling limit in sec.® and show that the "Russian doll" structure remains 
valid. The differential operators in that "scaled" nested Russian doll structure remain 
equivalent to the symmetric power of a singled out second order differential operator 
(corresponding to the modified Bessel function). In contrast, in the scaling limit, the 
direct sum of operators decomposition structure is lost, and we explain why. 

The unexpectedly simple expressions for the form factors fl^\ of sections 
([5]), and the corresponding remarkable differential structures, may be used to obtain 
many further results. We display some of these results in sec®. Recalling that, 
when A = 1, the Ising correlation functions C{N,N; 1) satisfy Fuchsian differential 
equations [19j with an order that grows with N, it is quite natural to inquire whether 
there are any other values of A for which C{N,N;X) will satisfy a Fuchsian linear 
differential equation. One such family of A is motivated by the work of Cecotti and 
Vafa |25| on TV = 2 supersymmetric field theories where they encountered A extensions 
of the Ising correlations in the scaling limit [21] with (m and n are integers) 

A = cos(7rm/ri) (18) 

^ The A-extensions ^ and JTOj are analytic at t ^ in t for T < Tc and, when T > Tc, analytic 
in t for N even, and in t^/^ for A'^ odd. 

J For the equivalence of linear differential operators, see |22l 1231 [24] . 



Holonomy of the Ising model form factors 



6 



Indeed, we have found that for n = 3, • • • , 20, the functions C{N^ N; A) satisfy 
Fuchsian linear differential equations whose orders, in contrast with those of the A = 1 
equations [TF], do not depend on N. More importantly, we find that these solutions 
are actually algebraic functions of t, associated with modular curves. 

We conclude, in sec. ([7]), with a discussion about the significance of our results on 
the factorization of multiple dimensional integrals. 



2. New integral representations for the fj^\^s 

The form factor expressions for C{M, iV) of [TOl [TTl [H [lH [M] are obtained by 
expanding the exponentials in and fS]), in the form given in [9] as multiple integrals 
and integrating over half the variables. The form of the result depends on whether the 
even, or odd, variables of [9J are integrated out. For the general anisotropic lattice, 
one form of this result is given, for arbitrary M and iV, in [14] . When specialized to 
the isotropic case the result is 



IM,N - ^ ^» Jm,N 

where s denotes sinh(2ii'), and where 



(19) 



^ ' jU-. 2 7: j^, 211 Vli sinh7„. 

x( n h.K)' (J[ COs(n ^n) (20) 



with 



l<i<k<j 



1 /I \ 1/2 

Xn — S ^ — cos 0, 



((s +--cos0„)2 -l) , (21) 

/ 1 \ 1/2 

sinh7„ = ((«+-- cos0„)2 - ij , (22) 



1/2 

2(x,Xfc)i/2 sin((0, - (/)fc)/2) 



(23) 



(7) 

In this work, we obtain the expressions of /jy V not by setting Af = in the 
results of [9], but, rather, from the representations of C{M,N) as an A^-dimcnsional 
Toeplitz determinant with elements 

^m,n — ^m — n (24) 

27r Jo V(l -aie-*«)(l -a2e»»)/ 

with 

ai = and a2 = s^^ (25) 
for the diagonal correlation C{N,N), and 



ai 



(l + .Y/^-.).(ii±^!^) and 



a2 = 1 + s 



((l + „Y'^-„)^((i±i!p±i) (26) 
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for the row correlatior0 C(0, N). Our method is to foUow Wu's paper [Tj in the 
framework of the general theory of Toeplitz determinants. 

For T < Tc, let us first recall (3.15) of Wu's paper [3, which reduces, for the 
diagonal correlations C{N,N), tc0: 

{1-t)-^/^ ■CiN,N) - 

^ ((1-^2 0(1 -«2/e))"'^' (27) 



1 



(2 7r)2 

d^'-r^^- ((1-^2 0(1- «2/n) 



1/2 



Comparing with ([7]) we see that the second term in ([27l) is /^^jy — JV- 
Performing the change of variables ^ = zi and ^' = 1 /z2 , deforming the contour 
of integration for both zi and 22 (one has to consider only the discontinuity across 
the branch cuj^ running from to a2), and rescaling zi and Z2, in, respectively, 
Xi — zi/a2 and X2 = 22/0^2, we obtain : 

fNNit) = F^vUt) = [' x^x, f X^dx2 



.Ti(l-X2)(l -tX2) ^ 

2:2(1 - a:^i)(l -txi) 



{\~txiX2)''^ (28) 



Similarly, when T > Tc, the leading term for G^^^ is given by equation (2.29) 
of 0: 

— 1 /" z^^^ 

((1 -il/2^)(l _ tl/2^-1)) 

which, after deforming the contour of integration to the branch cut, and scaling 
z — t^l"^ X, becomes 

(■(1) u\ - r<(i) 



/Ar,Ar(0 — GjV,Af(^) 

tN/2 ,1 



X 



^-i/2(1-.t)-i/2(i _a;t)-i/2rf^ 







,A^/2 (1/2)^ , 1. 



= t'''^-^^-2F,[-,N+-;N + l:t) (30) 

where 2^1(0, b; c; z) is the hypergeometric function [30] . 

The fuU expressions for F^^"^ for T < T^, and ^.^^ and 0^"+^^ for T > T^, 

can be obtained by following the iterative procedure based on (2.9)-(2.16) of [S] to all 
orders, just as the full expressions for Pj^^j^ and G^j"^^"* with M 7^ of (2.9)-(2.16) 
of [9] are obtained, in sections 3 and 4 of |9], by following the procedure of Cheng 

J Although in this paper we take the particular values of ai and 02 corresponding to C{N,N), 
our results are also, mutatis mutandis, applicable to the correlations C{0,N) and to the triangular 
lattice with 01 and 02 given by 1281 . 

^ To be precise, note that Wu considered in his paper, the C(0, A^) correlations. From the definition 
II24I I of the entries in the Toeplitz determinant one can consider the diagonal correlations C'{N, N) 
with the replacement I I25I I instead of I I26II . 
§ For T < Tc, 02 = t^l'^ < 1. 
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and Wu j8j to all orders]]]]. The details of the "iterative procedure" will be presented 
elsewherqj. These are certainly implicit in the paper of Jimbo and Miwa [2D] , however, 
we have not been able to find a reference where they are explicitly written out. 

When the low temperature expansion of sec. 3 of Wu [7] is performed to all 
orders, we find that (94) holds with 



F. 



(2«) 



N,N 



(-1)"+! 1 Zfdz, 



n T^^^ (31) 



X 

3 



n 



(2^)» J f \ 1 - z,zj+i 

{l~ a2Z2j){l~ a2/z2j) ^ 



V(l -a2Z2j-l){l -OL2lz2j-x), 

from which, after deformation of integration contours and rescaling, one obtains, for 
T < Tc, the following new integral representation of F^j^^{t) : 

= ^ ^ L X (32) 

1 2n N 



T-r dXj ^ " / X2j-_l(l - X2j){l - tX2j) \ 

jJi 1 J-jL \x2j{l -X2j_l)(l -tX2j-l) 



Similarly for T > Tc the expansion of sec. 2 of Wu [7] is performed to all orders 

,(2n) 



and we find that ^ holds with F^^^ given by (|32p and 



^(2n+l) _ ^ iv» 1 



n( zi z2„+i n 1 - Zjzj+i 

n+l 

X Y[ ((1 -a^^Z2j-i)(l -a2"V^2j-i)) 



■n. 



1 /2 

n((l -«2"'^2,)(l -«2"V^2,)) (33) 



X 

Changing variables and deforming contours, we obtain : 

^JV(2n+l)/2+2n 



^1 2n+l 



X 



Jo fJl ' XiX2n+l f}^ 1 

n+l / \ 1/2 

a;2i-i 



1 



xn 



(1 - a;2j_i)(l - tx2j-i) 



II The full expressions for F^^^ and G^^^j^ can also be obtained by performing the procedure 

based on (2.9)— (2.16) of [9] to all orders just as the full expressions for F^"lf and G^"jy^' with 
M 7^ of (2.9)-(2.16) of [9] are obtained in sections 3 and 4 of 0, by "cycling" the procedure of 
Cheng and Wu [8] to all orders. 

I The first step in that calculation is to consider the ratio C+{N, N,t)/C— {N, N,t). This will be 
detailed elsewhere. 
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1 /2 

\[({\ -X2j){l ~tx2,)/x2j) (34) 



X 



The form factor expressions are then obtained by expanding the exponentials. 
Thus we find, for T < Tc, that the form factors in © read 



X 



n 



a;2j-i(l-a;2i)(l -ia;2jO 



1/2 



7=1 V2:2j(l-a;2i-i)(l -ix2j-i) 

X J]^ J]^ {1 -tX2j-lX2ky^ 
l<j<7i l<k<n 

X J]^ (a;2i-i - 2^2fe-i)^(2^2j - a::2fe)^ (35) 

l<j<k<n 

and, for T > Tc, the odd form factors in read 



„(2n+l). ^ _ ,((2n+l)7V/2+n(n+l)) ^ ^ 

„1 2n+l n+l 

X / n ^fe' n ((l-^2j)(l -ia;2j)a;2i 
fe=i j=i ^ 

ri+l 

Y\ ((1 - 2^2j-i)(l - *2:2j-i)a:2j-i 



X 

i=i 

^ n n (1 -^a;2j-ia:2fe)"^ (36) 

l<i<»+l l<k<n 

X J]^ (.T2J-1 - a;2fc-l)^ ]^ {X2]-X2k)'^ 

l<j<k<n+l l<j<k<n 

where the last product in ([36]) has to be taken to be equal to unity for n = 0, 1. We 
note that the factors l/(n!)^ and l/(n! (n+ 1)!) in ([SS]) . and ([36|) . arise because the 
integrands are symmetric functions of the variables X2j and a;2j-i, separately. This is 
to be contrasted with ([^0]) where there is no separation in the odd and even integrals 

In the simplest case the previous integral representation p6p gives f\^\{t) defined 
by pop where one recognizes the Euler representation of an hypergeometric function. 

Do note that the {G'^^^^^\ ^jv"i Ar_(_i) decomposition in ^ is not unique. In 
contrast, the form factor expressions (j35p . (|36p are unique and well-defined. 

It is tempting to try to "bridge" such new integral representations ([35| , ([36]) with 
integral formulas like pO|) . or other integral formulas one can find in [71[8|,i9^,il0 t ill t i29j. 
getting ([55]) . from these other integral formulas after some changes of variables, 
or from partial integrations on a subset of variables in order to reduce (4n) integral 
formulas into (2n) integrals. We have not been able to do this. Basically we have two 
kinds of drastically different formulas: the ones emerging from Fredholm determinant 
expansions that naturally yield integral formulas with integrands that are algebraic 
functions of the self-dual variable w = s/(l -I- s^)/2, and the ones emerging from 
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"isomonodromic" calculationtQ that naturally yield integral formulas with integrands 
that are algebraic functions of the modulus k of the elliptic function (or the variables 
s or t\ for T > T^, k = , t = s^) and break the duality s 1/s. Do note that, in 
the case of the isotropic lattice, the formulas in [21 [H] are naturally integral formulas 
with integrands that are algebraic functions of the self-dual variable w = s/ (l + s^)/2. 
It is only in the scaling limit that these integral formulas look like integral formulas 
with integrands that are algebraic functions of s, or fc, or t. The " ^-integral formulas" 
of the second kind ([SS]) . ([36]) naturally produce series expansions in the hypergeometric 
functions 2^1, while the "w-integral formulas" of the first kind naturally generate 
series expansions in the hypergeometric functions 4^3 : see for instance, all the series 
calculations we obtained in [151 [HI HZl IHl in the holonomic analysis of x^'^'' ^-^d 
X^'*-'. The we consider are particular and, consequently, can be written, for fixed 
parameters, in terms of 2F1 for which quadratic transformationfl^ take place |30j : 

2F1 (a, b; 2b; 16w^^ = (37) 

(1 + ti/2)2a . (^a,a~b+ 1/2; b + 1/2; t) 

We have not been able to prove equality of the two kinds of formulas, so, instead, 
we have resorted to comparison of their series expansions. We have performed series 
expansions of our new integral representations in the variable t (see secQ), and found 
that they agree with the expansion of ([20l) . In the next section, we will see that they 
also agree with the coefficients of powers of A, h2j{N, N){t) and h2j+i{N, N){t), in 
the series solutions of (ITT]) - p5|) . 

Our new integral representations provide a "closed enough" representation of the 
A coefficients of the various A extensions C{N, N; A), the form factors. We will use 
the simplicity of these new integral representations in the next sections. 



3. Series solution of the sigma form of Painleve VI 

In this section we study the series expansions of the diagonal correlations C{N, N) 
starting from (jlip . the sigma form of Painleve VI. From order by order series analysis 
of the solutions of the equation |[TT|) . we show, when N is integer, the existence of a 
one parameter A extension of C{N, N), that actually identifies with the previous A 
extension of a Toeplitz origin in [7] , or of a Fredholm origin in 9J . 

We begin by considering some remarkably simple solutions of (jll[) . which exist 
for all N (not necessarily integer). In particular, consider the A'^-dependent second 
order hypergeometric differential operator |19| : 

It has regular singularities at t — 0, t — 1 and t — 00 with respectively the critical 
exponents (±7V/2), (1/4,1/4) and (l/4±iV/2). Denote by h any solution of ([55)1 . 
and consider the T > Tc expression (IT4| for a{t) : 

, s , ^. dlnh 1 
a{t) = t{t-l)-^ -- (39) 

X For instance formulas similar to formulas (13) and (14) of |20l . 

§ For a = b = 1/2 this is the Landen transformation | 30| on the complete elliptic integral K. 
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The hypergeometric differential equation l|38p. when written in a{t) given by 
takes a "Riccati" for ni: 

16 i (i - 1) • cr' +16-0-2 - 8 (i - 1) • CT (40) 

~ {2N ~1){2N +l){t-lf ^ 

For generic N ^ it can be verified that <j{t)^ given by (j39[) . is actually a solution of the 
sigma form of Painleve (llip . where h is any linear combination of the two solutions 
of ([38|) which, for generic N, read 

/± = t±^/2.(l-t)i/4. 2^1(1/2, 1/2 ±A;';l±iV;<) (41) 

and, for integer N, are /+ and 

. (1 _ i)i/4 . 2Fi([l/2, N + 1/2], [1], 1 - t) (42) 

We recognize from (|30l) that : 

/+ = (l-t)^/'-7T^-/ff^ (43) 

When t ^ the leading behavior of is /+ ~ t^^"^ which, if we make the 
normalization 

has the required behavior (|15|) for the high-temperature two-point correlation function 
C+{N,N) as t ~ 0. When the series expansion of is compared with the series 
expansion of C+ {N, N) [19] , we find that 

C+iN,N) = (l-t)i/4./^^)^ (45) 

_|_ (l/2)jv ■ ((3/2)Ar)2 ^3JV/2+2 , 

16 • r(7v -(- 2) • r(7v + 3)2 ' 

Thus, the first N + 1 terms of C+{N,N), and (1 - tf/^ f^^]j^, coincide. The 

coefficient of t3-^/2+2 q^u be considered as an "initial condition" needed to complete 

the characterization of the high-temperature two-point correlation function C+ (iV, N) 

seen as a solution of the Painleve VI equation (fTT|) . 

For T <Tc one notes that 

, dln((l-i)i/4) t , , 

.^t(t-l).^i-^-- (46) 

is a trivial solution of pTjl and that 

C_(^,iV)= (l-t)V4 + M(^f /)!)" ■ ^"^^ +•■■ (47) 

where the coefficient in front of is the initial condition defining the two-point 

correlation function C- ( A^, in the low-temperature regime jlQl [3T] . 

The relations (j45|) and (|47|) strongly suggest that, in order to analyze solutions 
of (|lip. we should introduce the following form for the low temperature expansions 

oc 

(iV, N)= (1 - <)i/4 . (l + ^ cfc • t'^) (48) 

fc=i 

% Note that, 1401 1. the Riccati form of II38I I. is also covariant by the Kramers- Wannior duahty JTHJ. 
This is a quite surprising result for a differential equation associated to f^^^{N, N), a form factor 
one could think to be specific of the T > Tc regime. 
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and, similarly, for the high temperature expansions : 

oo 

C+(7V, N)^ (1 - tf' ■ /i,^)^ + (1 - t) . . ^ 4 . (49) 

fc=i 

These expansions are not the most general solutions of (fTT|) because we have required 
that the solutions have the correct behavior (fT3| and (fTSll at t =0. These forms yield 
a one parameter family of solutions. 

We consider first the low temperature regime and use the form in (fTT|) to 
determine the coefficients recursively, order by order. When N is not an integer 
this recursive procedure gives the unique solution Ck = for all k. Thus the solution 
(1 — ty/^ is the only solution of the form (|48p . 

However, when N is an integer, we find that, while = for < iV, the 
equation which generically would determine cjv+i is automatically satisfied for all 
values of cn+i- The coefficient cjv+i can be specified arbitrarily and provides the 
second "initial" condition needed to specify a unique solution of (jlip. 

For all k > N + 2 the order by order procedure uniquely determines Ck as a 
polynomial in term of the free parameter cn+i- More specifically, the term (cn+i)"" 
first appears in the coefficient Cn{N+n)- Recalling (|47|) we see that if one writes the 
free parameter cn+i as 

(l/2)jv ■ (3/2)jv 
4-((7V + l)!)2 
the order by order solution to pT|) reads 



C_(iV,iV;A) = (1-0'^'- (^1 +E ^'"•^2„(iV)j (51) 

where /i2„(iV) for t - 0. This A extension C^{N,N;X) reduces to 

the low-temperature two-point correlation function C-{N,N) when A = 1, and to 
(1 — i)i/4 for A = 0. Using (|5ip in pTjl we have obtained the low-temperature series 
expansions of h2n{N) for a large set of values of the integer N, and found that these 
series expansions actually agree with the series expansions of (2n)-multiple integrals 
defined in sec. [H 

A similar order by order expansion can be carried out for the high temperature 
case. The corresponding coefficients dk can be deduced recursively, order by order. 
When N is not an integer the recursive procedure gives the unique solution dk = 
for all k. For non integer N we see that ft,jv = (1 ^ t)^^* ■ f^^^ is the unique solution 
of (HI]) of the form dMl). 

However, similar to the case for T < Tc, wc find that, when N is an integer, the 
coefficients dk are equal to zero for k < iV + 1, and that the coefficient dN+2 is a free 
undetermined constant. The coefficients dk, for k > A'^-l-S, are polynomials in dN+2- 
The term (dAr+2)" first appears in the coefficient rfn(Af+n+i)- Thus, recalling (|45p . we 
see that if we set 

_ (1/2)^ ■ ((3/2)^)2 

'^^^^ - ^ ■ i6.r(iv + 2).r(iv + 3)2 

the iterative solution to pip may be written, in the high temperature regime, in the 
form 

oo 

C+{N,N;X)= {l-ty/*-Y,^^"-h2n+i{N) (53) 

n=0 
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where 

hi{N) = f^^]^ (54) 

and where h2n+i[N) ~ j-n(N +n+i) t ^ Q. This reduces to the high-temperature 
two-point correlation function C_|_ {N, N) when A — 1 and to the hypergeometric 
function (1 - • for A = 0. Using (pjl in (fTT|) we have obtained the high- 

temperature series expansions of /i2j+i(A^)'s for a large set of values of the integer N , 
and found that these series expansions agree with {2j + l)-multiple integrals f'j^''^^^ 
defined in sec®. 

We have also performed low, and high, series expansions for the {N, N) defined 
by equations (4.2) in [T3] (see also (PU)) ). and we also found that these series identify 
with the one of h2j{N) and h2j+i{N) with the normalization : 

h2,{N) = C^^{N,N), h2,+i{N) = (55) 

It would be most satisfying if these identities could be demonstrated analytically. 

All the previous results confirm that these various A extensions identify and 
actually verify pT|) . The (log-derivative) of the A extensions C{N, N;X) satisfy the 
same (sigma-form of ) Painleve VI equation pT|) as the original diagonal spin-spin 
correlation, the boundary condition dependence coming from the original diagonal 
spin-spin correlation boundary condition. Even if some might consider that this result 
is not mathematically proved, it is clearly an exact result of experimental mathematics, 
based on an accumulation of large computer formal calculations. 

4. Fuchsian linear differential equations for /^^^(i) 

In previous studies on the Ising susceptibility [T31 [ini [13 [TH], efficient programs 
were developed which, starting from large series expansions of a holonomic function, 
produce the linear ordinary differential equation (in this case Fuchsian) satisfied by 

the function. In order for these programs to be used to study the ^'s we need 
to efficiently produce large (up to several thousand terms) series expansions in t of 

the /^^jy's- We have done this by use of both the integral representations ([35l) . ([36]) 

(7) 

and the representations of /j^ jv terms of theta functions of the name of elliptic 
functions, presented in 14J. 

We obtain the Fuchsian linear differential equations satisfied by the (diagonal) 
form factors f j^ j^ for j < 9. The analysis of these linear differential operators shows 
a remarkable Russian-doll structure similar to the nesting of (the differential operators 
of) the x*^^^'s found in [151 [16l [ITl [18] . Specifically we find that the expressions fl^^, 
/w AT' /^JV' /wjv actually solutions of the hnear ODE for f^j^^, and that f^^^, 
AT' /a/ TV' /a/ AT ^"^^ actually solutions of the ODE for /^^^. In addition, we find 
that all the linear differential operators for the /jy ^'s have a direct sum decomposition 
in operators equivalent to symmetric powers of the differential operator corresponding 
to f^\- Consequently, all the /^'tv's can also be written as polynomials in terms of 
the complete elliptic integrals E and K. The remainder of this section is devoted to 
the presentation of these results. 
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4-1. Fuchsian linear differential equations for f^j^"^^"^ 

The linear differential operator Fq{N) which annihilates fl^\ has the following 
factorized form : 

Fg{N) = Lio{N) ■ LsiN) ■ Le{N) ■ U{N) ■ L^iN) (56) 
where the differential operators Lr{N) are of order r. The first two read: 

9 2 1 - 1 1 1 



U{N) = L4,o - ■ L4,2 + Ye TT (58) 



with: 



(2<-l) , 24H^-241t + 46 
L4 n = Df^ + 10 • Dt^ + ^ ^ ^ • Dt^ 

{t-l)t 2 (t-lft^ 



(2t- 1) (122^2 - i22t + 9) 81 (5<-l)(5i-4) 

(< - 1)^ t3 "^16 i3 (i_ 1)3 



5 Dt^ (23 - 32t) 9 8-17< 

^4,2- 2 — ' 2(^-1)^3-^* "8(t^Ty7J (^'^^ 

The expressions (or forms) of Lq{N), Ls{N) and Lio{N) are given in [Appendix A| 
The linear differential operators F2„+i(7V), which annihilate fj^^i^^'' for n = 0, • • • , 3, 
are such that: 

Fr{N) = Ls{N) ■ Le{N) ■ LiiN) ■ L2{N) 
F5iN) = U{N) ■ U{N) ■ L2{N) 

Fs{N) = U{N)-L2{N) (61) 
F^{N) = i2(iV) 

Thus we see that the differential operator for /^'J^f '^'^ rightdivides the differential 

operator for Z^'^^'' for n < 3. We conjecture that this property holds for all values of 
n. We thus have a "Russian-doU" (telescopic) structure of these successive differential 
operators. 

4-2. Fuchsian linear differential equations for f^j^^^ 

The linear differential operator Fg,{N) (corresponding to f^\) has the following 
factorized form: 

FsiN) = Lg{N) ■ Lj{N) ■ L^iN) ■ L^iN) ■ ii(7V) (62) 
where the linear differential operators Lr{N) are of order r. The first two read: 

Li{N) = Dt, (63) 
, (2t-l) , (2-15t+14t2) 

8i2-l5i + 5 /Df 1 



2 (i - 1)^ t2 V i2 ^3 
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The expressions (or forms) of the hnear differential operators L^{N)^ Lj^N) and 
Lg{N) are given in [Appendix A[ 

Similarly to (|6ip there is a Russian-doll (telescopic) structure of these successive 
linear differential operators : 

Fe{N) = Lj{N) ■ L^iN) ■ LsiN) ■ L,iN) 
F^iN) = L5{N) ■ LsiN) ■ Li{N) 
F2{N) = L3{N)-Li{N) 

Fo{N) = L,{N) (65) 

Again, we see that the linear differential operator for /^'^ rightdivides the linear 

differential operator for /^'jj for n < 4. We conjecture that this property holds for 
all values of n. 



4-3. Direct sum structure 

Not only do the linear differential operators Lj{N) have a factorized Russian doll 
structure, but we have found that they also have a direct sum decomposition when 
the integer N is fixed. To illustrate this direct sum decomposition, the corresponding 
linear differential operator for f^'^pf reads: 

F3{N)= Li{N)-L2{N) = Ah{N)®L2{N) (66) 

where L2{N) is the linear differential operator for f^^^j^ and the fourth order operator 
Mi{N) is displayed in [Appendix B| for successive values of N . One remarks on these 
successive expressions that the degree of each polynomial occurring in these linear 
differential operators M^^^N) grows linearly with TV. 

As a further example consider f^^\N^ N), where we find that the corresponding 
linear differential operator decomposes as: 

= Le{N) ■ L4{N) ■ L2{N) = Me{N) ® M^iN) ® L2{N) (67) 

where L2{N) is the differential operator for fj^^, M4(iV) is the previous fourth order 
differential operator, and the sixth order operator Mq{N) has again coefficients whose 
degrees grow with N for successive values of N. There is nothing specific to fj^ pf 

and f^\ ■ similar results hold for all the /^]v'^' being even or odd. 

In contrast with the Russian-doll way of writing the differential operators for 

(n) 

f)^ )y, the direct sum structure, as a consequence of this growing degree, cannot 
be written for generic N as operators with polynomials in front of the derivatives. 
This "non-closure" of the direct sum structure will have some consequences when 
performing the scaling limit of these differential operators (see sec.® below). 



4-.4- Equivalence of various Lj{N)'s and Mj{N)'s linear differential operators 
We find that the symmetric square of L2{N): 

. (2t~l) , (l-7t + 7i2) 
Sym\L2{N)) = Dt' + 3 ^-^^—^^ ■ Dt' + ^-j—^ ■ Dt 



^ The symmetric j-th power of a second order linear differential operator having two solutions /i 
and /2 is 
solutions. 



and /2 is the linear differential operator of order j + 1, which has , /j , . . . , /j 
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1 l-2t Ar2 ^ A/-2 

- 7^ 0— Dt - 7 — ^ 68 

and the linear differential operator L-i{N) are equivalent : 

L:i{N)-U{N) = V{N)-Sym^{L2{N)) (69) 
with the foUowing intertwinncrs : 

U{N) = {t-l)t-Dt^ +{3t~l)-Dt +1 + • iV^ (70) 

V{N) = {t-\)t-Dt^ + (11 t-5)-Dt 

(5t-l)(5^-4) (t-1) 2 

(i-l)i t ^ ' 

Similarly, with the symmetric cube of L2{N), we have the equivalence 

L4{N)-A{N) = B{N)-Sym^{L2{N)) (72) 

with : 

AiN)= it-l)t.Dt' +^(2t-l).I)t2 + l£^^liH±^.Di 

9 2f-l 9(t-l)-N'^ ^ 9(2t-l) 

+ o 7 7^ T ^ Dt -- ^ • 73 

8 (f - 1) i 4 i 8 i2 ^ ^ 

B(Ar) = {t-l)t-Dt'^ +y (2i-l)-£>t2 

21 (6-29t + 29t2) 9 (2i - 1) (125*^ - l25t + 16) 

^ T (t-l)i ^ 8 (f - if ^2 

9 (^- 1) ..2 ^ 9 (lOi- 9) 

- - -5^ '- -N^ -Dt - - -^^ ^ — '- ■ (74) 

4 i 8 t2 y I 

More generally, all the Lm{Nys are (m — l)-symmetric-power of L2{N). As 
a consequence their solutions are (m — l)-homogeneous polynomial of the two 
hypergeometric solutions of L2{N). 

Similarly, for the linear differential operators occurring in the direct sum, one 
easily verifies, for every integer TV, that, for instance, the M4(iV)'s are equivalent to 
the cubic-symmetric-power of L2{N): 

Mi{N)-Q{N) = S{N)-Sym^{L2{N)) (75) 

where, for iV = 0, 1, 2 : 

Q(0) = {t-l)t-Dt +t - (76) 

Q(l) = 2 {t-lf ■ Dt^ +2, {Z-7t + At'^){t-l) t- Dt"^ 

^ ^„ 9 41 9\ ^ 3 2i2-2i-|-l 
l2t^-2%e + -t--yDt +- , (77) 

Q(2) = \ {t-lf {3 + 8t + 3t'^) t-Dt^ 

+ ^ {l5-t-35t'^ + 15t^ +6t'^) {t-Vj-Dt^ 

1 (18 - 12 - 97 + 577 1^ - 738 i -|- 252) 
~ 24 t 

1 12 -I- 14 - 260 -I- 497*2- 314 1 + 24 
"16 P 



Dt 



(78) 
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As a further example, one can verify, for every value of the integer N , that the 
sixth order operator Mg {N) is equivalent to the fifth symmetric power of L2 {N) . The 
solutions of the linear differential operators M^^N) are also (to — l)-homogeneous 
polynomials of the two hypergeometric solutions of L2{N). 

As a consequence of this direct sum decomposition, the solutions /("^(A^, N) 
are (non- homogeneous) polynomials of the two hypergeometric solutions of L2{N) or, 
equivalently, /^^jy see dH])) and its first derivative. The second order linear 

differential operator L2 [N) is equivalent [12] to the second order differential operator 
Le 

Le = At- Dt'^ +4:Dt - (79) 
corresponding to the complete elliptic integral E. As a consequence of the 

(n) 

previously described direct sum decomposition, the /j^jy'^ '^^^ written as 

polynomial expressions of the complete elliptic integral E and its first derivative E' , 
or alternatively, E and the complete elliptic integral K . 

(2) 

Let us just give here a set of miscellaneous examples. For ff^ pf^ one has: 
2 tfo = {K -E)-K 
2 f^^l = 1 -3KE -{t-2)-K^ 
6 t ■ f2-2 = &t - {2 + Qt^ - lit) ■ K"^ 

- {15t- A) ■ K E -2{l + t)-E'^ (80) 
90iVi^3 = 135t^ ~ {I37t^ ~242t^ + 52t + 8) ■ K'^ 

+ (8 - 319 + 112 i + 16) • /sTi; - 4 (1 + t) (2 + 13 1 + 2) E^ 
3150 • fi^l = 6300 

- (32t^ + 2552i2 + 128 + 6440 1^- 1119U^ +464t) • 
+ (128 + 5648 - 14519 1^ + 1056 i + 576 1^ + 256) • E K 
-8 (1 + t) (I6t^ + 58t3 + 333i2 +58t+16) ■ E^ 
where E and K are given by (I17p . Other examples are given in [Appendix C[ 

Remark: All these remarkable structures are not restricted to diagonal two points 
correlation functions. Actually one can calculate various j-particle contributions 
/m^at off-diagonal two point correlation functions, and verify, again, that they 

are also polynomial expressions of the complete elliptic integrals E and K. For 
instance for T > Tc : 

C(2)(0,l)= I -l{l + .s')K (81) 

-IeK -i(.2-3)(l + .2).x2 

where s = sinh(2if). Other miscellaneous examples of such off-diagonal j-particle 
contributions are displayed in [Appendix D[ 

4-5. The elliptic representation of Painleve VI 

The results we have underlined in this section, namely the unexpectedly simple and 
remarkable polynomial expressions for the form factors /jv jv correspond to the 
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fact that the associated hnear differential operators are direct sums of operators 
equivalent to symmetric powers of the second order differential operator Le- We 
already encountered this central key role played by the linear differential operator 
Le, or the hypergeometric second order differential operator (|38p . in our previous 
holonomic analysis of the two-point correlation functions of the Ising model [19] . In 
order to understand the key role played by Le-, or equivalently operator ([38]) . it is 
worth recalling (see |32j . or for a review [33]) the so-called ^^elliptic representation^^ of 
Painleve VI. This elliptic representation of Painleve VI amounts to seeing Painleve VI 
as a "deformation" (see equation (33) in [33]) of the hypergeometric linear differential 
equation associated with the linear differential operator : 

C= {l-t)t- Dt^ +{l-2t)- Dt (82) 

One easily verifies that this linear differential operator is actually equivalent (in the 
sense of the equivalence of differential operators) with Le, or equivalently ([55]) . 
This deep relation between elliptic curves and Painleve VI explains the occurrence 
of Painleve VI on the Ising model, and on other lattice Yang-Baxter integrable models 
which are canonically parametrized in term of elliptic functions (like the eight- vertex 
Baxter model, the RSOS models, see for instance [M]). We will see, in sec®, another 
example of this deep connection between the transcendent solutions of Painleve VI 
and the theory of elliptic functions, modular curves and quasi-modular functions. 



5. The scaling of /^'jy 

(n) 

The scaling of the f}^ j^'s amounts, on the functions, and on the corresponding 
differential operators, to taking the limit N ~* oo and t — + 1, keeping the limit 
X = N ■ (1 — t) finite, or in other words, to performing the change of variables 
t — 1 — x/N, keeping only the leading term in N. Performing these straightforward 

calculations, the linear differential operators in t for the /^j^'^ where N was a 
parameter, become linear differential operators in the only scaling variable x. 

Calling _Fj"="' the scaling limit of the operator Fj (N) we find for j even that 

j^scal 
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(83) 



where 



^sca/ ^ 2 x^Dx^ + 10 x'^Dx^ - 2 (7 + 5 x^) Dx^ 

+ 2 (-16 + 13 x^) x^Dx^ + 2 {5-12x^ + 4:x'^) X Dx 
-10 + 8x2-24x'', 

= 2 x^Dx^ + 8 x'^Dx^ 2 {x ~ 1) (x + 1) x Dx ~ 2, 
Lf"' = Dx (84) 
and Lq'^'^\ are given in [Appendix E[ 
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Similarly, for j odd, we have 
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(85) 



L|™' = IGx^Dx^ + 96 ZJa:^ + 40 (2 - a;^) Dx^ 

+ 8 {x^ ~2) xDx +9x^ -8x^ + 16, 
Lf-^' = 4a;Da;2 +41)3; - x (86) 



^2 

and -C/fo"', ' are given in [Appendix E| 

" 3 



Thus, we see that the scaled operators F^'^"-^ have a "Russian doll" structure 



inherited from the lattice operators Fj (N) . 

Consider the linear differential operator corresponding to the modified Bessel 
function Bessel{n,x /2) for n = 0, namely 

B = Dx^ + — ~- (87) 
X 4 

We recognize, in this linear differential operator, the exact identification with the 
scaled differential operator F^'^"-^ = L^°-^ . We find that the symmetric square of the 
linear differential operator _B, and the scaled operator if' are equivalent : 

• {xDx^ +2Dx -x) = (88) 

{2x^Dx^ + 12x^Dx ~2x^ + 8x^)- Sym^(B) 

Similarly, the symmetric third power of the linear differential operator B, and the 
scaled operator Lf^' are equivalent, and, more generally, the symmetric j-th power 

^ scat 



of ([57)1 and the scaled operator L^™' are equivalent 



Lff, c Sym^B) (89) 
Recall that the differential operators Fj{N), corresponding to the form factors 
/i^ Af' '^^^ written as direct sums only when the integer N is fixed. At the scaling 
limit, this feature dissappears in the scaled differential operators Fj'^"''' which have no 
direct sums. Therefore while the scaling limit preserves the Russian-doll (telescopic) 
structure (see (I5T|) . ([55)) ) and also preserves the fact that the various operators in this 
Russian-doll (telescopic) structure are equivalent to symmetric powers of an operator 
(j87p which replaces the operator Le, the direct sum structure is lost. As a consequence 
the scaling of the /jv cannot be seen as simple polynomials of modified Bessel 
functions. 

(2) 

There is one exception that concerns fj^ j^- Its scaled linear differential operator 
pscai ^ YiQ,^ the non shared property of being equivalent to the direct sum of Dx with 
the symmetric square of (|87p . namely: 

pscal ^ j^scal^pscal _ Dx®SyTn^{B) (90) 

From this equivalence, one immediately deduces the expression of the scaling of the 
/jv w quadratic expression of the modified Bessel functions of x/2 which actually 
identifies with formula (2.31b)-(3.151) in [9]. 
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The occurrence of niodified Bessel functions, emerging from a confluence of two 
singularities of the complete elliptic integrals E and or from the hypergeometric 
function 2^11 should not be considered as a surprise if one recalls the following 
limit of the hypergeometric function 2F1 yielding confluent hypergeometric functions 
iFi. These confluent hypergeometric functions, iFi, are nothing but modified Bessel 
functions [5U] : 

2; 

iFi(a, 6; z) 2-^1(0, p,b; —) when: p ^ 00 

V 

F^r^|^-^^(-+i2.+l;2z) (91) 

Remark: It was shown, in sec.(j4]), as a consequence of the decomposition of 
their differential operators in direct sums of operators equivalent to symmetric powers 
of Lb, that the functions /^]v polynomial expressions of E and K functions. 
Therefore their singularities are only the three re^wZar points t = 0, t = \ and t = 00. 
The scaling limit {t = I — x/N , t 1, N 00) corresponds to the confluence of the 
two regular points t — and t — 00, yielding the, now, irregular point x = 00. The 
occurrence of irregular points with their Stokes phenomenon, and, especially, the loss 
of a remarkable direct sum structure, shows that the scaling limit is a quite non-trivial 
limit. 

Contrary to the common wisdom, the scaling limit does not correspond to more 
"fundamental" symmetries and structures (more universal ...) : this limit actually 
destroys the remarkable structures and symmetries of the lattice model^ 



6. Algebraic solutions of PVI for A ~ cos(7rm/ri) and modular curves 

The function C{N,N;X) is such that its log-derivative is actually a solution of the 

sigma form of Painleve VI : it is a transcendent function "par excellence" . However, 

a) 

the unexpectedly simple expressions for these form factors strongly suggest to 

try to resum the infinite sums (O, and ^TU\i . of form factors Z^''^, corresponding to the 
function C{N, N; A), and see if these transcendent functions could be "less complex" 
than one can imagine at first sight, at least for a set of "singled-out" values of A. For 
instance, are there any values of A 7^ 1 which share, with A = 1, the property that 
C{N,N;X) satisfies a Fuchsian linear differential equation ? 

Actually, introducing, instead of the modulus k of elliptic functions (for T > T^, 
k — s^), or the s and t variables, the name of the elliptic functions (see relations 
(5.7)-(5.11) in [H]), we have been able to perform such a resummation, getting, 
for arbitrary A, nice closed expressions for the C{N, N; A) for the first values of 
N, {N — 0, 1 , 2 • • •), as sums of ratios of theta functions (and their derivatives), 
corresponding to Eisenstein series, or quasi-modular forms. These results will be 
displayed in forthcoming publications. The simplest example corresponds to iV = 
where C-.{N, N; A) is just the ratio of two Jacobi 6*3 functions : 

C_(0,0;A) = ^4^' where: A = cos(u) (92) 

"3(0, q) 

§ These kind of results should not be a surprise for the people working on integrable lattice models, 
or on Painleve equations |35l 136) . 
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All these results strongly suggest to focus on u — 7rm/n (m and n integers) yielding 
for the possible choice of "singled-out" values of A : 

A = cos{nm/n) (93) 

Actually these special values (1^5)) of A already occurred in a study of = 2 
supersymnietric field theories [25 in a similar series construction of solutions of the 
Painleve V (or Painleve III for a ratio of functions) equation for the scaling limit of 
the Ising model [9]. 

We have begun to investigate this situation. When n — 3, • • • , 20 (and all 
the possible values of m, but a set of first successive values of A^), we have found 
that C±{N, N; X) do indeed satisfy Fuchsian linear differential equations but, unlike 
the equations found in [IQJ for A = 1, the order of the Fuchsian linear differential 
equations depends only on n and not on N. 

As examples of these Fuchsian linear differential equations, we found, for instance, 
that C_(Af, A^; cos(7r/4)), for A^ = 0, 1, 2, are annihilated, respectively, by 

= it-irt-Dt' +^-it-mt^2)-Dt +1 

4^'^ = {t-lft- Dt^ +{t- l)(5t - 2) . - 1^ + ^ 
4^/''l = (t - 8)(t - if t-Dt^ +{t- l){t^ -2t + 16) • Dt 



209 1^ 25t 1 



(94) 



256 16 2 
and that C_ (0, 0; cos(7r/3)) is annihilated by : 

4/3] ^ _ ^^3 ^3 . ^^4 -l){t- If e ■ Dt"" 

+ ^ (43^2 - 43t + 4) (t - 1) i • Dt^ (95) 



(2t - l)(247t^ - 247t - 80) • Dt 



35 



1458 ^ ' • 486 

These linear differential operators are of a quite different nature from the one depicted 
in sec.(|4|) which can be decomposed in direct sums of (operators equivalent to) 
symmetric powers of L^;. In contrast with the direct sum decomposition we have 
underlined previously, these linear differential operators are irreducible. However we 
do expect from sec. (|4.5l) a connection with elliptic curves. Actually, instead of a 
connection through the second order differential operator Le, or the hypergeometric 
second order linear differential operator ()38p . we have an even more striking link 
with the theory of elliptic curves. These solutions C{N,N; A) are actually algebraic 
solutions of Painleve VI, associated with modular cMrwe|§|. We found for n < 8 these 
singled-out Fuchsian linear differential equations, corresponding to algebraic solutions 
of Painleve VI, and beyond, directly these modular curves for larger values of n for 
which we do not have the Fuchsian linear differential equations yet. 

We first obtained these modular curves as polynomial relations P{a, t) = 0, 
between a and t, and we then found, in a second step, the polynomial relations 
P(r, t) = 0, between r = C±{N,N;cos{TTm/n)) and t. For instance, one finds that 
r = C_(0, 0; cos(7r/3)) is solution of a genus one algebraic curve : 

16t^^ -16t^ -8{t^l)t-T^ + t-{l-t) = (96) 

§ The occurrence of modular curves is pretty clear for N = from Il92|l . from the analysis of its 
invariance group, subgroup of the modular group. 
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or that T — C^{N,N;cos{n/'i)) is solution of genus three algebraic curve, for 
instance, for = 0, 

16t^^ +16{t-l)-T^ +t'^ -(t-l) = (97) 
the corresponding solutions being quite simple algebraic expressions : 

C_(0,0;cos(7r/4)) = 2-1/^(1 - + (1 - t)i/2]i/4 (gg) 

C_(l,l;cos(7r/4)) = 2-3/'*(l - + (1 - i)i/2]3/4 ^gg^ 

C_(2,2;cos(7r/4)) = (100) 

= 2-5/4(1 - t)l/16[l + (1 - t)l/2]5/4[5 _ _ ^)l/2]/4 

We give in Table 1, when available, the order of the Fuchsian linear differential 
equation for A = cos(7r/n), the degree and genus of the corresponding algebraic 
curve P(C_(0, 0;A), t) ~ 0, and the degree and genus of the algebraic cr-curve 
P(cr(0,0; A), t) = 0. 



n 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


14 


16 


18 


20 


ODE order 


4 


2 


12 


4 


24 


8 


















r-degree 


12 


16 


60 


48 


168 


128 




240 














r-genus 


1 


3 


13 


13 




41 


















(T-degree 


4 


2 


12 


4 


24 


8 


36 


12 


60 


16 


24 


32 


36 


48 


(T-genus 








1 





4 







1 




1 











Table 1. Order of the linear ODE, as well as degree and genus of the 
corresponding modular curve in r = C±(0,0;A) for A = cos(iT/n), when 
available. The corresponding degree and genus of the modular curve in (t(0, 0; A), 
when available. 



We found the following results on the polynomial relations P(r, t) = 0, between 
T = C-{N^ N] cos(tt m/n)) and t. These polynomials are actually polynomials of the 
variable p = r" for n odd and of the variable p = " for n even. This property is 
related to the invariance of the variable p under a subgroup of the modular group0. 
Let us denote (5(p, t) = the polynomial relation between p and t. We also found 
that the degree of the polynomial Q in p actually identifies with the degree in a of 
the polynomial P((t, t) = 0. Thus, the r-degree in Table 1 can be seen to be the 
(T-degree multiplied by n for n odd, and by 2 n for n even. The order of the Fuchsian 
linear differential equations for C_ {N, N; A) identifies with that degree in a. We 
finally found that the genus of the modular curve P(cr, t) = identifies with the 
genus of the t" (resp. T^")-modular curve Q{p, t) — : the genus corresponding to 
C_(0, 0; cos(7r/3))^, C_(0, 0; cos(7r/5))5 are respectively 0, 1 , the genus corresponding 
to C_(0,0;cos(7r/4))8, C_(0, 0; cos(7r/6))i2^ C_(0, 0; cos(7r/8))i6 are but the 
genus for C_(0, 0; cos(7r/10))^'' is 1. In contrast the genus corresponding to 
C_(0,0;cos(7r/6))^ C_(0, 0; cos(7r/8))8 are 1 , and the genus for C_(0, 0; cos(7r/10))i" 
is 5. 

For = 0, and only in this case, a large set of these algebraic curves (for 
instance (j96p or the modular curve for n = 7 in the previous table) are invariant 

^ See in particular Barth and Michel 1371 for further details on the jfoo(n, 2) modular curves and 
the characterization of the genus of modular curves from subgroups of 51/(2, Z). We will study 
C± [N, TV; A) from this modular subgroup point of view elsewhere. 
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under the t ^ \ — t symmetry : 

(t, cr, o-', cr") -> {l-t, -cr-1/4, cr', -cr") (101) 

This remarkable symmetry is, in fact, inherited from the covariance by (jlOip of the 
sigma form (|lip when iV = 0. 

A large set of algebraic solutions of Painleve VI (and associated modular curves) 
have been obtained by many authors [551 [5^ HUl BT[ [i^ . However, most of these 
results on algebraic solutions are for the canonical forniiJ of Painleve VI in terms of 
the variable y: 

d^y I ( I 1 I \ f dyV (l 1 I \ dy 



dt"^ 2 \y 1-y y~tj\dtj \t t-1 y - 1 J dt 
y{y~l){y-t) ( t-1 , , t(t-l) ^ 



t2(t-l)2 V ' ' y^ ' ' {y-iy {y-ty 

There are several sets of a, /3, 7, 5 which lead to the same equation [IJ for a. For 
T < Tc one such selQ of parameters of p02p , corresponding to the A^-dependent sigma 
form ((TT|) . is : 

<.^i(iV + i)\ p^-i[N-i)'. ..i. (103) 

It is interesting to make the connection between our results and those previously 
known algebraic solutions. Such a "dictionary" will be performed elsewher^, let us 
just give one simple example. The variable y being a rational expression of a and its 
derivatives (see [H]), the algebraic solution ((55)) with (|103p . becomes y = 1 — \J\ — t 
which is the well-known solutior{^ y = \/t (see [40]) under the change t ^ 1 ^ 
y ^ 1 — y, P ^ —7 and 7 — /3 which is a symmetry of (jl02p . 



7. Conclusion 



The diagonal Ising two-point correlation functions can be expressed (see for 
instance [HI [45]) as homogeneous polynomials of complete elliptic integral E and 
K. These diagonal Ising correlations are A = 1 subcase of their A-extensions 
C{N, N; A) we considered in this paper. By ([7]) and ([8]) these polynomials of 
E and K are also expressed as infinite sums of the form factors f]y/j~f's which, 
themselves, are polynomials of E and K. This yields a double infinity (M, N) of 
remarkable identities on the complete elliptic integrals E and K. Similarly, with 
the previous algebraic solutions for A = cos(7rm/n), one sees that an algebraic 
expression C{N, N;cos{TTm/n)) (associated with a modular curve) can be written 
as an infinite sum of polynomials in E and K. Each of these modular curves will 
provide a remarkable identity on the complete elliptic integrals E and K . 

Recalling relations like (5.7)-(5.11) of [14J, all these identities can also be written 
in terms of the nome of the elliptic functions occurring in the Ising model. These 
identities, now, become remarkable identities on some infinite Gaussian sums, or on 

II For N = this equation has been solved in terms of theta functions 1401 1411 142| . has dihedral 
symmetry and has a countable number of algebraic solutions, 
fj To be considered when comparing with |40|. 

^ We found that the SL{2, Z) subgroup for t" for n odd, (resp. t^" for n even) identify with the 
one for y (see [37] [40] . 

§ The solution y = y/i solves II102II for the parameters [a, /3, 7, 5) = (a, —a, 1/2 — 5, S). 
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series expansions of theta functions, or, for large enough values of iV, on Eisenstein 
series and other quasi-modular forms. We will describe, and analyze these identities 
in a forthcoming publication. 

The calculations displayed in this paper can be seen as successful explicit examples 
of factorization of multiple integrals, providing examples of explicit calculations of 
the new mantra that "nested sums are Hopf algebras and thus multiple Feynman- 
like integrals must factorize in terms of polynomial expressions of one-dimensional 
integrals" . For our j-particle contributions of the diagonal correlation functions, 
the C^-'-*(M, iV)'s, the fact that they are polynomial expressions of singled-out one- 
dimensional integrals (the complete elliptic integrals E and K) is understood in 
terms of direct sums of linear differential operators equivalent to symmetric powers 
of a singled-out linear differential operator. In the scaling limit, this direct sum 
structure, yielding polynomial expressions (that is the so-called "factorization of 
multiple integrals"), is lost: what remains is a Russian-doll structure of differential 
operators equivalent to symmetric powers of a singled-out differential operator. 

The problem of the factorization of multiple integrals is, obviously, an important 
one for Feynman-like integrals. It also occurs on various calculations of correlation 
functions of integrable models, like the Heisenberg spin chain, where multiple integrals 
also occur. These factorizations are obtained by Boos and Korepin [46l|47] by adding 
to the integrand a successive set of anti-symmetric integrands (these anti-symmetric 
integrands being chosen in such a way that their multiple integral is zero). The 
combination of the initial integral with these new integrands yielding expressions 
depending on less variables, thus reducing the n-multiple integrals to a (n— l)-multiple 
integral. More recently, Boos et al. [48] also deduced factorization of multiple integrals 
representing the density matrix of the Heisenberg spin chain: the key ingredient, in 
the emergence of such factorization, is a functional identity on the integrand, this 
relation coming from the Bethe ansatz integrability of the model. The factorization 
of some multiple integrals can probably be seen as a consequence of some "Yang- 
Baxter integrability" , it seems, however, to occur beyond this narrow framework. The 
Feynman-like integrals, where such factorization of some multiple integrals occurs, 
are not arbitrary holonomic expressions. What are the (more or less integrable) 
constraints one must impose on holonomic integrands such that their multiple integrals 
exhibit factorization, remains a fascinating open question [49j . A key point we have 
tried to promote here is that, instead of trying to calculate multiple integrals where 
the integrands have no free parameters, that is to say that the multiple integrals are 
just constants [50| . we perform calculations on multiple integrals where the integrands 
do depend on one, or many, parameters. We can then use the holonomic structure. 

In short: it is simpler to get multiple integrals that depend on one variable than 
obtaining their evaluation at a given value on this variable. This is typically a Yang- 
Baxter view point : it is easier to solve an integrable model with a spectral parameter 
that enables to describe the Yang-Baxter structure, than trying to solve that model 
for a given value of that parameter (quantum groups, knot theory, etc.). It is easier to 
solve the anisotropic Ising model than the isotropic one, and, similarly, it is easier to 
consider multiple integrals that depends on a variable, than evaluating constants |49j 
(polynomial expressions of C(3)), C(5)j •■•) corresponding to these multiple integrals at 
a given value of that parameter: this way of looking at the problem enables to see the 
emergence of highly non trivial algebraic structures on linear differential operators, 
that are a very efficient and powerful tool of experimental mathematics, and other 
formal calculations, to study factorizations of multiple integrals. 
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Appendix A. Differential operators Lj{N) 

The linear differential operators Lj{N) have the following form: 

n=0 k=0 y'^y^ ^)) 

where uq = (j — l)/2 for j odd, and no = j/2 for j even. The polynomials P^l{t) 
are of degree fc in t. 

Appendix A.l. P!^l(t) 

P^^^ = 1, p^^l = 40 1 - 20, Pp' 2^ = -563 1 + 558 1"^ + 118, 

„(5) 4291 10169 , „„„„ , „„„ 

P^V = 1 1^ + 3320 1^ - 220, 

0,3 2 2 

P^^l = 80 + 10848 - 16978 1^ + 8180 - 2227 1, 

P^'^^ = 4 (85 - 1139t + 3672*2 _ 4259^3 _^ leoot'') t, 

P^(_5q) = -5, P^^l = -9H + 59, Pi^a* = -469*^ + 626t - 181, 

P^^^ = 144- 840 i+ 1368*2 - 656 ^^ P^^^ = 4, P^^l = 16t- 16 



Appendix A. 2. Pi^}{t) 



<J=1, Po^?= 70^-35, Pg^ = j^t^-r^t + 413, 

Po^^3 = 2 (2 i - 1) (5912 1^ - 5912 1 + 979) , 

„(6) 2410523 4 2410523 . 3200163 o 98705 „„„ , 



py°^ = —(2t-l) (3585925 1^ - 7171850 1^ + 4326453 1'^ - 740528 1 + 19600) , 
P^^^ = ^t{t-l) (48841 1'^ - 97682 1^ + 63549 1'^ - 14708 1 + 784) , 
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(6) _ 34799 2 43231 6133 

„f6l 88609 , o 96849 16691 

^3 = —t"" + 41823 e - t + 

P^(*^] = [t - 1) (94091 - 146523*2 + 67548 f - 9216) , 
p(6) _ 259 (6) _ 1917 3159 

Pi,«)=^(407i- 272)(t-l), P^o^=-^ 

Appendix A. 3. Pi]l{t) 

P^J^ = 1, Pq(_7 = -56 + 112 Po^5 = 5012 - 5026 1 + 1148, 
P^^l = -10736 + 797274- 174373 + 115544^^ 
Pq(J2 = 46172 - 5487364+20429534^ - 29752444^ + 14728284^, 
P^^l = -78640 + 1605642 4 - 9634279 4^ + 23975501 4^ 

- 26144958 4* + 10305440 4^ 

P^J = 29160- 16160784+^1^^4^ 

136608085 , 511207495 , ..o....^ r 
^ 4^ + 4^ - 111249042 4^ + 36334360 

P^J^ = ^4^59940 - 16650374+ 11865715 4^ - 36308026 4^ + 544662944^ 

- 39393900 4^ + 10951200 4^) , 

P^(JJ = -14, P^(^ = -966 1 + 57 A, 
p[^l = -24712 + 29686 4 - 8248, 
p[^l = -290812 4^ + 530547 4^ - 299013 4 + 51188, 

p[^l = -1561136 4'' + 3851903 4^- 3309480 4^ + 1156221 4 - 136440, 

r,(7) .^r.r.r.r. 22166415 , 18989235 4 11893977 , 

P\'>^ = 129600 4^ + 4^ + 4^ 

^'^ 2 2 2 

2902725 5 
4- 30281044^, 



^ppendia; P„°fe(4) 



2 

Pj^J = 49, Pi_7 = 16864 - 1254, P^^l = 178874^ - 270264+9679, 
P2^3^ = -22761 - 133569 4^ + 577534^ + 982534, 
Pg^Jg) = -36, Pj^^ = 324 - 3244 

(8), 



P^^^ = 1, Pg(®) = 168 4 - 84, P^^^ = 11697 4^ - 116974 + 2730, 
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P^^l = 2 (2 i - 1) (109862 1^ - 109862 1 + 21881) , 

^(s) ^ 77675835 _ 77675835 ^3 ^ 108450015 _ 7693545 ^ ^ 3g^3^^^ 

Pq(_S5) = i (2t - 1) (257365313t^ - 514730626 + 340542345 

- 83177032 i+ 6033464), 

r,(8) 135579123 4351723053 , 31150612733 , 
Po^e' = 2610671 1 + 1^ 1^ 

47738959467 . 55357772589 . 15912986489 . 

r H r H r, 

16 16 16 

P^^^ = i (2t - 1) (16309728941 1^ - 48929186823 + 54824769942 
' 8 \ 

- 28100895179 1^ + 6440184015 1^ - 544600896 1 + 8016008^ , 
P^^l = t{t-l) (719580625 f - 2158741875 1^ + 2496751275 1'^ 

- 1395599425 1^ + 383051976 1^ - 45042576 1 + 1308736^ , 

P[^l= -21, P[^l = -2352 i+ 1365, 

„(8) 414555 2 483843 
-Pi 2 = : — t + — : — * - 33315, 



4 4 

P[^l = -2290461 1^ + 4034358 1^ - 2237787 1 + 386664, 

„(8) 426526863 , 504203159 , 845513895 , 36865265 

P\ l = 1^ + t^ t^ H \ 

16 8 16 2 



- 223043L 

61()r)8()181 - 7342474719 , 4139827129 . 

- 4 * + 16 * 8 * 

4378085671 „ 

+ — 1^ - 67155042 1 + 6072033, 

16 

p'f^ = -1^ [t - 1) (449304249*5 - 1168884874i^ + 1134316077 

- 509448428 1^ + 105774112 1 - 829440o) , 

(8) _ 987 (8) _ 15993 _ 22299 

-^2,0 — g ' 2,1 ^2 4 ' 

(8) _ 2933043 2 4128099 696405 
^2.2 - 16 * 16~*+~8~' 
(8) _ 7002915 3 _ 14949545 2 10236397 _ 4465707 

2.3 - i * i * + i * § ' 

pfl = _ 1) (2179797*3 - 4103797*2 + 2457908* - 468864) , 

' 128 

(8)_ 3229 (8)_ 21963 76827 

^3,0 - > ^3,1 - 4 16 ' 

^iS=-^(6607t- 5032)(*-l), = ^ 
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Appendix A. 5. P^^^it) 

Pq^ = 1, Po!i = -120 + 240 1, P^^^ = 5796 - 24546 1 + 24516 i^, 

P^^^ = -145528 + 991701 1 - 2099751 1^ + 1396208 ^^ 

P^^l = 2045004 - 20325858 1 + 69369177 1^ - 97902648 1^ + 48749364 i*, 

P^^^ = -16074560 + 225525578 1- 1125696965 + 2565535675*^ 
- 2714936962 1^ + 1079617840 ^^ 

(9) 37765468163 , 63124281313 , 

Po,6 = 66126712 - 1333788966 1+ t^ 1^ 



+ ^16166^483 ^4 _ 45734526046 1' + 15125870712 1^ 

r,(9) o oor> ^r>o..„ 78510959875 , 374049548401 , 
= -118102672 + 3823928460 i H 1^ 

- 471178501099t'^ + 646530989251^"^ - 455734056216*^ 

7 fg) 138902716891 , 
+ 128906004992 *^Po\s^ = 47071232-4139526516*+ 1^ 



, 3495148889889 . 
484196478836 1^ + - 3539969007392 1 



5 



2 

+ 4054878125399 1^ - 2448333931344 1'^ + 604418968592 1^, 
Po^9 = 80*^9561344 -427020633* + 4937178194*^ - 26308505171*^ 

+ 76760779797** - 130255661861 *^ + 128108854250 *'' 

- 67626000000 f + 14796800000 **) , P^^^ = -30, 

P^^l = -5082* + 2898, P[^^ = -352662 *2 + 404466* - 110238, 

P^^^ = -12963996*^ + 22438245*^ - 12306435*+ 2123604, 

p[^l = -271930980*^ + 631696597*^- 523169724*^ + 181823257* 

- 22193940, P^i = -3245449704*^+1^555^^*4 

21116262613 3 11093266991 n 2731643299 

*^ + * * + 125146416, 

2 2 2 

p[^^ = -20342103432**^ + 71951600804*^- 100847772344*4 
+ 71245445309 *^ - 26573340926 1'^ + 4930067225 * - 354631488, 

P^^^ = 406425600- 363525018400*^- 6762200560* + 46013156464*2 
+ 214239244800 *^ - 51475353600 *^ + 324098542224 *" 

- 162977694704*3, P^^^ = 273, P^^} = 29490* - 19650, 

P^^^ = 1217265*2 - 1636902*+ 528465, 
P^^i = 23917695 1^ - 48731759 + 31834675 * - 6637935, 
P^^l = 222934641 - 612371540 + 607781638 - 257369288 * 
+ 39119361, P0 = 734599360 * - 2367055040 1"^ - 2730725376 ** 
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- 87745536 + 785703936 + 3664705792 1^, 
P^^^ = -820, P^^l = -46428 1 + 37212, 
P^^^ = -839284*2 + 1358312* - 535156, 
P^^^ = 2455552 - 9463296* + 11831808 - 4814848 ^^ 
Pio^ = 576, P^} = -9216 + 9216* 

Appendix A.6. P^]^^^) 

1, =-165 + 330*, Hf= 11286 -l^*+l^*^ 

P^^P = 440 (-1 + 2 *) (4400 *2 - 4400 * + 947) , 

r.(io) oK^^^^^o 2289151821 „ 1601770797 , 

Py?' = 9053979 - 85922628 * + 1^ *^ 

>J,4 8 4 

1601770797 , 

P^^°'^ = i (2* - 1) (^27291921049 *'^ - 54583842098 *3 

+ 37651935321*2 - 10360014272*+ 946138408^, 
p(-)= 907059937- ^^^^^^^^^^^+^^^^^^^^^^^^2 

10567258749853 , 17677263640199 . 14872361118327 



^3 ^ -- ^'^ _ 



32 32 32 



4957453706109 



g 

+ 32 ^ ' 

P^^°^ = i (2*- 1) (9240801571631 *•' 

' 8 V 

- 27722404714893*^ + 32221157315067** - 18238306771979*^ 
+ 5175369000414*2 - 676616400240*+ 30201789392^, 

P^f = 7010881775 - t + ^^^^^^'^^^^'^^^ f 

2 16 

320491542697265 . 16973016403001045 , 

*3 H ** 

16 256 

8195527196507945 , 18373024724608855 » 
*^ H *'' 

G4 128 

5532574254401525 ,7 5532574254401525 § 



t + ^rr. i > 



64 256 

p(i°) = (2*- 1) f2925753951778285*^- 11703015807113140*^ 
256 ^ ^ V 

+ 19103849088522126*^- 16350991940670388*^ 

+ 7838575034697949*" - 2079015276577248*^ + 280113055050736 *2 

- 15268105688320*+ 14855363712o), 

-Pq.io = (8079810760125 *^ - 32319243040500 *^ 
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+ 53723369995078*^ - 48052759343484 + 24975072368117*^ 
- 7567996044344*3+1257903576048*^ - 96158271040* 

1833995520^, P[]o'> = = -10032* + 5643, 



+ — J, - i,u - ^ , 

(10) _ _ 4116057 2 4646169 _ 629409 

A,2 - I t + - * — 



„(io) 116107101 . ^^r..r.r.r.r. 1 107291349 18577449 

P| o ' = + 98646966 * + , 

2 2 2 

(10) _ 62944154655 4 71575322887 3 117218126643 2 

32 16 32 

10187763521 630658425 

+ 1 ' 

(10) _ 163658617341 5 1868493139019 4 1024185567025 3 
- \ * + * § * 

1072700805259 2 .^^„^„.o.o^ 3132363327 
+ *^ - 16676053189*+ ^ , 

(10) _ _ 32591037777225 g 112104944104795 ^ _ 154297814907493 4 

n,6 - 64 * + 64 * 64 * 

108234764940653 . 20292699152369 o 

H r r 

64 32 
958545911705 35460987675 

+ 1 * 4 ' 

(10) _ 110727690476325 7 446374095368415 g 

^'^ 32 32 

740746002842197 . 163197508324913 , 82052435490193 , 
^ ** *3 

32 8 8 

93651745559635 . 1748299590545 105491089125 
+ 32 ' 4 4 ' 



32 ' 4 4 ' 

' =-T^ (*-l)fl530421397125*^- 5554611547375*^ 
1024 ^ ' \ 

8232623167111 f"^ - 6428537243541 f!^ + 2844044623496 
71 1552088080 f + 93740238400 * - 5138022400^ , 



- 71 1552088080 f + 93740238400 * - 51380224 

(10) _ 4389 (10) _ 364353 _ 469491 

- 8 ' ^2,1 - 4 i 8 , 

(10) _ 194340135 ^ 251951799 19747365 
2.2 - ^2 * 32 8 ' 

„(io) 830796045 , o 1027545163 208821765 
Pi = *3 - 406585301 1^ + * , 

^.•^4 4 4 

(10) _ 490847729943 4 645152215543 3 1232916142207 2 

128 64 128 

63246182437 4689902523 

— — * + — 5 — ' 



10 o 
(10) _ 2314309478331 5 15325677165945 4 1231151157175 3 

64 128 8 

_ 12252529589789 2 230003188957 _ 26662508757 

128 * + 8 * 8 ' 

-^2,6°^ = ^ (*-l) (95338644413*5-286574346250*^ 
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+ 332882516705 1^ - 186665173556 1^ + 50609686768 1 - 5337817088^ , 

(10) _ _86405 (10) _ _ 1071807 3230007 

3.0 - 32 ' 3,1 - 4 16 ' 

(10) _ 616031665 ^ 934953233 173220155 

- 64 * ^ 64 ^ 32 ' 

„(io) 4731457901 , 2715664857 , 8133623529 1981312349 

' = t^ H t^ t H , 

32 8 32 32 ' 



P^''^> = -—— {t-l) ( 191281007 399820191 + 271766508 1 
- 60099968V 



^(10) _ 1057221 (10) 23053617 . 41642109 



256 ' " 128 256 ' 

= lii (892447^- 727072) (t - 1) , = 



Appendix B. Direct sum structure 

We display the fourth order differential operator M^^N) introduced in sec. 14.31 for 
successive values of N: 

, (2t- 1) (2t^ - 2i + 3) o 
MAO) = Dt^ +2--^ -^-Dt^ 

^ ' {t^ -t + l) {t-\) t 

1 (-73<+14+102t2-58i3 + 29i4) , 



2 {t'i - t + I) {t - ly 

1 (2 1 - 1) (5 - 10 13 + 27 f2 - 22 i + 2) 

2 (V - i+ 1) (t- l)3t3 
1 -2t^ + 'i2t^ -'ilt + 'i 



Dt 



16 (t2-t+l) (i- 1)3^3 



M4(l) = i^t^ +2- Z?f3 +1 £2 2 

^' {t-l)t-P4 2^t~lft^Pi 



+ $ +4. ^° 



where 



{t ~ 1)3 ^3 . 16 _ 1)3 ^4 . 

Po = 256 1^ - 560 + 312 - 143 1^ + 227 - 72 i - 16, 

Pi = 64 - 856 1^ + 2826 - 4087 1^ + 2978 1^ - 1098 + 182 i - 8, 

P2 = 208 + 7807t2- 14253 13 + 12412 i*- 4624 + 448 ^^ 

P3 = 64 - 556 1'' + 1225 1^ ~ 1078 1^ + 396 1 - 48, 

P4 = 16 + 209 1^ - 120 1^ + 16t^ - 120 1, 



M4(2) = 



1 



P3 

t (t - 1) P4 



P>t3 



Pi 



2 t3 (i _ 1)3 



P't + — 



1 

2 

Po 



P2 



1)' P4 



■Dt^ 



16 i4 _ lY 
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where the corresponding P,;'s read : 

Po = -1344 + 10752 + 139321 1^- 721147 + 1888781 

- 3452437 1^ + 4219535 1^ - 3184189 1^ + 1330028 1^ 

- 202384 1^ - 34048 1 + 7168, 

^ 448 + 4256 + 56658 i^- 519911 + 1502563 1^- 2077796 
+ 1426525 1^ - 372047 1^ - 39536 1^ + 5418 1^ + 14336 1 - 896, 
P2 = 4928 37632 1394407 + 4810853 1**- 8001289 
+ 6880493 + 415793 - 2881207 + 16128 + 11648 
+ 174818 

F3 1344 1^ + 1568 1^ + 65828 1'^ - 382102 1^ + 760238 1^ - 702181 

+ 302183 1^ - 46627 + 1568 1 - 1792, 
P4 ^ 448 1^ + 448 + 16513 f - 81242 1^ + 127675 1"* - 81242 1^ 

+ 16513 + 448^ + 448 

Appendix C. The form factors /jv jv 

In order to check all the results displayed in this paper, we have performed a large 
number of series expansions. Even the series expansions obtained recursively, order 
by order, from the sigma form of Painleve VI pT|) in sec. 3, were checked against 
series expansions obtained independently. Some were based on extremely large series 
expansions, not in s or i, but in the nome of elliptic functions (see (5.7)-(5.11) of fT4^), 
others were obtained from series expansions with hypergeometric functions coefficients. 

Actually our new simple integral representations (I35I1 . (f36|) are of a great help 
to produce large series expansions for the quantities fj^"^{t) and fj^^j^^Ht). This 
amounts to expanding only the (1 — ^2:2^-1X2^)"^ term in (|36p . Recalling the Euler 
representation of the hypergeometric functions [30j : 

1 

x''-^ {1 - xy-''-^ {1 -xtydx 



one can rewrite, alternatively, these integral representations (|35p. (|36p expansions of 
^^"i fNlv^\'^)^ nested sums of products of hypergeometric functions. By 
expanding the factor (1 — txiX2)~'^ in a power series in t we obtain 

F(-l/2, N + j + 1/2, N + j + 2;t) F(l/2, N + j+ 3/2, N + j + 2;t) 

The series expansions for h2j{N, N){t)^s and the h2j+i{N, N){t)'s agree with the 
series expansions for the C^{N, A^)'s and with the series expansions for the /^^^'s 

and /^jv^^'s- In this Appendix, we display the /^''^r's and /^'^v^^'^ some j and 
some N . 
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Appendix C.l. f^^^^j^ and f^^^j^ 

The fl^^^s are given in the text by (|¥T|l and (|44p. The /^''^'s are given exphcitly as 
a function of K and E in the text. 

Appendix C.2. f^^j^j 

The /^^^jv's read for N = 0, • • • , 4: 

6 /^J^o' = K-{t-2) -3K^E 

6i^^^ /f.i = 4:{K- E) -6K^E - (2i-3) +3KE^ 

18 1 /fa^ = 7 + 2) - 14 + 1) + 24 

+ 3 (2 _ 11 f + 2) if^ „ 3 (-^2 _ 2) ^3 ^ 3g (-^ _ ^) 

270 ^^/2 f^^^ = -30 (8t2 + 7t + 8) t-E + 30 {At^ + 3t + 8) K 

~ (72 - 158 1^ + 189 1^ - 156 1 + 8) iC^ 

+ 6 (24t'^-108t^ + 29i2-6t + 4) i;^^^ 

+ 3 (232f'^-lll<2-i80t-8) S^Js: 

+ 4 (t + 1) (2 + 103 1 + 2) tE^ 
M2h0t^ ffl ^ 975 (3t + 4) (8^2 _ 5t + 12) t^i^ 

- 7800i^(t+ 1) (6i^ -< + 6) 

- (^16216 - 32109 1^ + 4218 + 38472 1^ - 38064 + 3264 1 + 128^ 

+ 3 (l0832 f - 43424 1^ + 4925 + 13248 1^ - 10112 + 3328 1 + 128^ E 

- 48 (4 1^ - 2885 + 939 + 1510 + 1792 + 212 1 + 8) 

+ 16 (8 + 216 + 4893 1^ + 5464 + 4893 1^ + 216t + 8) E^ 

Appendix C.S. f'^^ 

Some of the f'^\'s read: 

24 /^^J = 4 (iC - i;) • - {2t-3)K^ -6 K^E + 3 isT^i;^ 
2Af[^} = 9 -30i^i; -10(i-2)if2 

+ (t^ - 6 1 + 6) if * + 15 X^i;^ + 10 (t - 2) if^ E 
72 1 ■ f!fl ^ 72t-32 (l + t) E^ 

-16 {2 + 6t^ -lit) -16 {15t- 4) KE 

+ (24t3-98i2 + 113i-36) 

+ 12 {9 + t) E^ K +3{71t- 60) K'^E^ 

+ 2 (66 + 74 1^ - 157 1) EK^ - 24 E^ 

1080 1^ f^^l = 22 (8 1^ - 319 + 112 1 + 16) KE 
-88 (l + t) {2t^ + 13t + 2) E^ 
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+ (957 - 3646 1^ + 4230 1^ - 1488 1 - 8) i^^ 

+ 8 (46 1^ + 51 + 543 1 - llO) E^K 

+ 3 (16 - 72 + 2537*^ - 2704t + 272) K'^E^ 

- 6 (8 - 903 1^ + 1934 1^ - 988 1 + 40) EK^ 

+ 24 (13 + 13*2 - 28 1) E'^ 

- 22 (137 1^ - 242 1^ + h2t + ^) K"^ + 2025 1^ 

Appendix C.4- /^jv 
We give some f]^\'s: 

120 1^/^ /|;\' = 64 • - £;) 

- 20 (2 1 - 3) - 120 EK^ + 60 E^ K 

+ {4:t -5) {2 t-:i)K'' + 15 (2t-3) X'' £; +A5K^E'^ -IbK^E^ 
360 * /^J^2^ = 149 (t + 2) K ~ 298 (t + 1) ^ 
+ 720E^-90 (t^-2) ii"^ 
+ 90 (2t2_iif + 2) EK^ + 1080 {t-l)E'^K 

+ (5*^ + 28*2 -90 1 + 60) - 10 (t^ - 16*2 + 24* -4) EK^ 

- 5 (32 *2 - 179 i + 122) K^E'^ - 30 (19 i - 29) K'^E^ - 360 E'^ K 
5400t^/Vi;3 = 792 (4i2 + 3i + 8) i-Jf - 792t(8i2 + 7i + 8) 

-40 (72i*-158i^ + 189t2-l56t + 8)iir^ 

+ 240 (24i*- 108*^ + 29*2-6* + 4) EK'^ 

+ 120 (232 - 111 *2 _ 180 * - 8) ^ 

+ 160 (t+ 1) (2*2 + 103* + 2) E"^ 

+ 5 (96 f"^ - 520 + 1310 - 1589 + 800 * - 88) 

+ 5 (424 - 2488 * + 5051 12 - 4962 + 2008 - 192 1^) E 

- 5 (1984 - 9228 + 9423 12 3272 1 + 816) E^ 
+ 5 (784 - 4104 i + 1 1697 12 - 6056 1^) E^ 

- 40 (2 *3 + 738 *2 _ 567 * + 47) E^K 
+ 360 (*2 - 28*+l)£;^ 

Appendix C.5. f^\ 

Some /^^jv read: 

720 = 225 - 259 (* - 2) ii:2 _ j^j k E 

- 105 Jf^i;^ + 525 ii'2i;2 + 350 (i - 2) E +35 {6 + t'^ -6t)K^ 
-21 {6 + t^ -6t) EK^ -105 {t- 2) E^K^ 

-(i-2) (t2-lOi+10) 
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2160^/^^2^= 2160i - 544 (15i-4) 

- 1088 (1 + t) -E"^ - 544 (et^- lK + 2) -X^ 

+ 50 (24t3-98<2 + 113f-36) - 1200£;* +600 (9 + K 
+ 150 (71 1 - 60) K'^E^ + 100 (66 + 74 - 157 1) S 
+ ?,QOKE^ - 15 (235i- 264) if^£;^ 
+ 3 (720 - 1889 1 + 1490 1^ - 344 1^) ■ E 

- 90 (21 + t) E'^K^ - 45 (92 + 74 - 173 1) • E^K'^ 

- 3 (32 1^ - 220 1^ + 504 1^ - 467 ^ + 150) 

All these fl^\ displayed when expanded have their leading coefRcients starting 

as given in ([551) ^'^'^ Let us give some as a series to show the magnitude of 

the numerical coefficients involved. The series expansion of fl^\ for the first values 
of N reads: 

= 



37 1 



10 



N = 1 
= 2 
A^ = 3 

Af = 4 



1073741824 8589934592 
7ii2 21ii3 



4398046511104 2199023255552 
21 1^^ 19215 



1125899906842624 144115188075855872 
10395 84315 



18446744073709551616 18446744073709551616 
2335905 t^i 166783617 



75557863725914323419136 604462909807314587353088 



Appendix C.6. fYl, f^l and /{\^ 
for the 

5040 /f 1^ = 2304 -{K- E) 



Here we give the /^\r for the other values of j = 7, 8, 9 and A^ 1: 



- 784 {2t-3)-K^ -AYOAK^E + 2352 AT i;^ 
+ 840 {2t-3) K* E - 840 E^ 

+ 56 (4t-5) {2t-3)K^ + 2520K^E^ 

-28 (4t-5)(2i-3) EK^ -210 (2t-3) E^ 

- (32 - 156 + 228 1 - 105) K'^ + 105 AT^i^* - 420 K^E^ 

645120 = 11025 - 38748A:£; - 12916 (t - 2) AT^ 

+ 29610 K^E^ + 1974 {t^ -6t + 6) K* + 19740 {t - 2) A'^ E 

- 8820 E^ - 84 (< - 2) (i^ - 10 1 + 10) 

- 1764 (i^ - 6 1 + 6) A:^ A - 8820 (t - 2) K* A^ 
+ 945 /sT^'A'* + (t^ - 20 + 48 - 56 1 + 28) AT^ 

+ 36 (< - 2) {t^ - 10 < + 10) A:^ A + 378 (t^ - 6 i + 6) AT^A^ 
+ 1260 (t - 2) E^ 
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362880 j(9) = 147456 ■ {K - E) 

+ 157U0 KE^ - 52480 {2t-3)-K^ - 3U880 K^E 

- 65520 K"^ E^ + 65520 {2t-3) K'^E 

+ 4368 (4 1 - 5) (2 1 - 3) fC^ + 196560 K^E"^ 

+ 12600 - 50400 K^E^ ~ 3360 {4t - 5) {2t - 3) E 

+ 120 (105 - 228 1 + 156 - 32 1^) K'^ 

- 25200 {2t-3)K'^E^ + 630 {U - 5) i2t - 3) ■ K'' E^ 

- 945 + 4725 i^^ + 3150 (2 1 - 3) E^ 
+ 45 (32 1^ - 156 + 228 i - 105) E 

+ ( 128 - 960 1^ + 2460 - 2572 1 + 945) 

Appendix D. Miscellaneous off diagonal j-particle contributions 

We display here some off diagonal j-particle contributions. 

2s2 • C(2)(0,2) = 2s2 -2{l + s^)s^-K 

+ {2s^ + -2) ■ -{s-2){s + 2) KE -2E'^ 

8s4-C(2)(0,3) = s2 (8 + 27^2 + 8/) 

-24s2(l + s2) -E - 2 8^{\^s^) (-4+13s2 + 8s'') K 
+ (1 + 5^) (8s8 + 7s<5 + 3s4 -8s2-8)-i^2 

+ 4 (4 + 6s2 + 7s^ + 6s'') - 8 (s'' + l) • 

18s^-C(^)(0,4) = 36 (2 + 8s2 + 13s'' + 8s^ + 2s^)s2 

- 24 s2 (1 + s^) (6 + 18 + 11 - 12 §2 _ 8) ii- 
-48s2(i + s2) (7^4 + 15^2 + 7) E 

+ (s^ + 2) (72 + 144 - 60 - 200 5^ + 5^^ + 62 - 16) 
+ (64 - 408 - 576 s'^ + 591 + 1088 + 336 3^°) £; K 

- 4 (8 - 93 - 200 s'^ - 93 s*^ + 8 /) E^ 
8s2.C(2)(l,2) = 7,s2 -4 (l + s2) -i? 

- 2 (1 + 5^) (.2 - 2) ■ , 
+ (l + ,s2) (4 + 3s2-5s^) 

+ 4 (s^-3s2-3) Eii' +8^2 
6s4-C(2)(l,3) = 95^ -4 (l + s^) (54 + 3^2^1^.^ 
+ 4 (l + s^) (1 +3,s2-,s4) - if 
+ (10 + 8 _ 2 - 8 s^' - 5 s^) • if 2 
+ (-24 - 32 - 13 s'' + 16 + 4 s^) • Ji- 
+ 2 (7 5* + 12s2 + 7) 
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Appendix E. Differential operators in the scaling limit 

The differential operators Lg'^"', • • • , introduced in sec. Oread: 



Ll^''^ ^ GAx^Dx^ + 320x^Dx^ - 16 (48 + 35 x^) Dx'^ 
+ 32 x^ {91 ~ 4) Dx^ +Ax^ (848 - 1788 x^ + 259 x^) Dx'^ 

- 20 a; (80 - 380 x^ + 383 x'^) Dx 

- 225 x^ - 2480 x^ + 1600 + 17580 x'^ 

^""^ = 4 Dx'^ - 56 (3 + x"^) x^ Dx^ + 8 (41 + 84 x^) x* Dx'^ 
+ 4 (69 - 810 + 49 x'^) x^ Dx^ 

- 8 (251 - 971 x'^ + 372 x'^) x^ Dx^ 

+ 4 (275 - 2116 + 4212 x^ - 36 x^) x Dx 

- 1100 + 2832 - 35280 + 1152 x^ 

L|™' = 256x*i:»x^ - 1024x^W - 5376 (2 + x^) x^ Dx^ 
+ 256 (334 + 399x2) x^ Dx^ 

- 32 (4040 + 26304 x^ - 987 x") x^ Dx^ 

- 64 (4216 - 57384 x^ + 12027 x^) x^ Dx^ 

+ 16 (76688 - 537424 x'^ + 482478 - 3229 x^) x^ Dx^ 

- 16 (48400 - 598032 x^ + 2328262 x^ - 60013 x^) x Dx 

+ 774400 - 3342592 x^ + 72498272 x" - 4879248 x^ + 11025 x^, 

Ll"""^ = 4x^L»x^ + 480x^i:)x^ -24(5x^-961) x^Dx^ 
+ 8 (71315 - 1092 x^) x^ Dx^ 
+ 12 (645013 - 19478 x^ + 91 x'^) x^ L'x''^ 
+ 8 (6985303 + 4920 x'' - 354291 x^) x'^ Dx^ 
+ 4 (44460417 - 3774790x2 + 108828 x'' - 820 x^) x^ Dx^ 

- 4 (443021 + 5872124 x^ - 382676 x'^ + 9216 x^) x^ Dx^ 

+ 4 (576 x^ - 16128 x^ + 94812 x"* + 9265148 x^ - 268975475) x Dx 
+ 1024 (36 x^ - 2019 x'^ + 66804x2 _ i254400) 

Lj;^"' = 1024 x^° Dx^" + 168960 x^L'x^ + 8448 (1368- 5x2) x^ Dx^ 
+ 11264 (37880 - 399x2) ^7 jj^7 

+ 4224 (^2194904 - 44164 x2 + 133 x"*) x^ Dx^ 

+ 128 (946138408 - 30108276 x2 + 259215 x"*) x^ Dx^ 

+ 32 (29025917984 - 1305848840 x2 + 21377796 x^ - 86405 x^) x" Dx"^ 

+ 64 (60403578784 - 3569603544 x2 + 92712956 x^ - 1057221 x^) x^ Dx^ 

+ 4 (^1794785734400 - 134201812672 x2 + 5056843872 x'^ 

- 110074968 x*^ + 1057221 x^) x2 i:>x2 
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+ 972 (3056659200- 330174912x2 + 18778592 - 694968 a;*^ 
+ 18375 -xDx - 893025 
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